Two-fluid behavior of the Kondo lattice in the 1/N slave boson approach. 
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It has been recently shown by Nakatsuji, Pines, and Fisk [S. Nakatsuji, D. Pines, and Z. Fisk, Phys. 
Rev. Lett. 92, 016401 (2004)] from the phenomenological analysis of experiments in Cei-^La^ Coins 
and Celrlns that thermodynamic and transport properties of Kondo lattices below coherence tem- 
perature can be very successfully described in terms of a two-fluid model, with Kondo impurity and 
heavy electron Fermi liquid contributions. We analyze thermodynamic properties of Kondo lattices 
using 1/N slave boson treatment of the periodic Anderson model and show that these two contribu- 
tions indeed arise below the coherence temperature. We find that the Kondo impurity contribution 
to thermodynamics corresponds to thermal excitations into the fiat portion of the energy spectrum. 

PACS numbers: 71.27.+a, 75.30.Mb, 75.20.Hr, 72.15.Qm 
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I. INTRODUCTION 

Anomalous behavior of heavy fermion (HF) metals has 
attracted enormous theoretical and experimental interest 
over the past 30 years since their discovery Due to the 
almost localized nature of f-electron bands, electrons in 
these materials have strong Coulomb interactions, and 
exhibit a variety of very unusual properties, including 
heavy effective mass of carriers (with m ~ 100 — 1000m e ), 
non- Fermi liquid ground state properties |2j, nontrivial 
magnetic order (proposed for UPt3, URu 2 Si 2 [jl, 0, 0)j 
and unconventional superconductivity, observed in most, 
if not all of HF superconductors (see, for example, Refs 
0, for a review) . Recent theoretical and experimen- 
tal interest in these materials has been focused on the 
quantum phase transitions. Quantum critical anomalous 
power law behavior with temperature of thermodynamic, 
magnetic, and transport properties^ is observed in the 
vicinity of a quantum phase transition in HF compounds 
such as CeCu 2 Si 2 |l, CeCu 5 qAu .i 0, YbRh 2 Si 2 |l, or 
CeCoIns and Celrlns [13, [H|l3. The surprising experi- 
mental discovery is that the quantum criticality in heavy 
fermion metals is also unconventional, i.e., it does not 
follow the usual Hertz-Millis|H [l| scaling analysis of 
a magnetic quantum critical point. Instead, the fluctu- 
ations remain local, as it was demonstrated in the in- 
elastic neutron scattering experiments Q, which st rong ly 
suggests a Kondo-type local quantum criticality0, Il5| - 

The most interesting question remains the formation 
and the origin of heavy mass in these materials. Nakat- 
suji, Pines, and Fisk 12](NPF) have recently shown that 
our understanding of Kondo lattices is incomplete. NPF 
propose that formation of heavy mass in heavy fermion 
materials (in particular, in Cei- X L& X Coins and Celrlns), 
which happens at a crossover temperature T* , can be de- 
scribed by a phenomenological theory similar in spirit to 
Landau theory of superfluidity [lal . The NPF two-fluid 
phenomenological model^, Il7j states that thermody- 



namic, magnetic, and transport properties of Kondo lat- 
tices can be described as a sum of two independent con- 
tributions, one involving the HF liquid, the other a lattice 
of weakly interacting Kondo impurity centers. The rela- 
tive fraction / of the coherent HF state plays the role of 
the order parameter, which develops below the coherence 
temperature, T*. 

We demonstrate below that the same type of phe- 
nomenology directly follows from the lattice analogue of 
the slave boson broken symmetry ansatz of Read and 
Newns(18j and Coleman 19]. The order parameter, the 
slave boson average, gives rise to the coherent heavy 
fermion liquid component, while the contribution from 
almost flat f-electron bands accounts for the Kondo impu- 
rity component. Although no actual symmetry breaking 
occurs in this case, an analog of Landau two-fluid model 
in superconductors arises in the slave boson language. 
The heavy Fermion liquid contribution corresponds to 
the bose-condensate of slave bosons, while the Kondo im- 
purity contribution is a complete analogue of the normal 
part [Ig . 

While the slave boson language provides a very attrac- 
tive physical explanation of the NPF phenomenology [l2j . 
in does not completely agree with it. We provide a de- 
tailed comparison of the slave boson results with the NPF 
phenomenology in the summary. 



II. THE SLAVE BOSON METHOD 

Theoretical approach to understanding the physics of 
HF usually starts with the Kondo lattice, or more general 
periodic Anderson p(i| model Hamiltonian, 



H — Hn + H\. 



(1) 



where 
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^ fl\fi\,fii f il- 



ka 



(2) 



Here the creation and annihilation operators for the f- 
electrons, f\ a and carry the site index i, and there is a 
Coulomb interaction at each site for the f-electrons. The 
operators c ka and c ka correspond to delocalized Bloch 
states. The hybridization term Hy accounts for the c — 
/ hybridization between the f-electrons and the Bloch 
states: 



constraint 24] for every /-site i: 



E 

i.k.a 



T t ikRi n\ i T7* —ik-Rr t r 

V k e 'fi a c ka + V k e *c' ka U 



(3) 



Both Anderson and Kondo models have been solved 
exactly for a single f-electron site (for a review, see 
Ref.[2ll Even though the Bethe ansatz solution 

is by no means easy, most of our basic theoretical knowl- 
edge about Kondo systems is inferred from it. Unfortu- 
nately, exact solution cannot be obtained for a periodic 
lattice, where correlations involve several conduction and 
f-elcctron spins. Nevertheless, it can be treated quite 
successfully using various approximate methods, such as 
the non-crossing approximation, the 1/N approach |2 1| . 
or dynamical mean field theory [23(. 

The large-TV slave boson approach was developed in 
the mid-80-s by a number of authors[H [H H ill for 
single impurity Kondo and Anderson models, and was 
later applied to lattices [H HHH |2^. The result of 
this approach for the single impurity models has been 
shown to approximate well exact results at both low and 
high temperatures [3(j. Following these references, we can 
rewrite the periodic Anderson model Eq.Q in the slave 
boson language. Since the on-site Coulomb repulsion U 
is very large for f-electron materials, it is usually taken 
as infinite. Then the creation and annihilation operators 
for /-electrons need to be applied with special projection 
operators, which project out doubly occupied /-electron 
statesj24fl . A convenient way to deal with this Hamilto- 
nian |l8t H^. l25j is to introduce a new slave boson field 
bf , which creates a hole on site i, and to rewrite the 
Anderson Hamiltonian in a way which allows the 1/N 
expansion in the number of orbitals. As a result, the 
infinitc-C/ Anderson model takes the following form: 



H = 



E [ ekC km C krn + ^o/fc m /fc« 
km 

( C lmf<irnbl_ k + &.C.J 



(4) 



k.q.m 



Here m is a spin index, which runs from 1 to N (i.e., we 
assume N degenerate f-levels and N conduction bands), 
and the sums for k and q run over the Brillouin zone. 
The absence of doubly occupied states is ensured by a 



Qi=Y.fimfim + b\b i = l. 



(5) 



The large- N slave boson formulation and the 1 /N expan- 
sion around it have many well-known difficulties. Nev- 
ertheless, the results of this model capture the low tem- 
perature heavy Fermi liquid properties of Kondo lattices 
rather well, and agree with many experiments |2l| . The 
mean field approach starts with a replacement of the slave 
boson operators with their averages [T^. I24I l28j , 



(h) = 



(6) 



After taking into account the constraint Eq.JSJ), the re- 
sulting Hamiltonian then can be written as: 



H = E \_ e kC km C km + tffhmfkr, 



+ (Va)(ft m c km + h.c.) + (e f - E )(a 2 - 1)] , (7) 

where e/ is the renormalized /-energy level. The parame- 
ters e/ and a are fixed by the mean field equations, which 
are obtained by minimizing the free energy. 

A serious well-known limitation of the 1/N slave bo- 
son approach arises from the necessity to satisfy the con- 
straint Eq.JSJ). According to quantum mechanics, the 
phase of the boson field is conjugate to the number oper- 
ator for it. Thus, a state in which the boson has a definite 
phase, such as one introduced by the usual slave boson 
mean field ansatz Eq.JSJ), results in infinite uncertainty 
for the boson number operator, and necessarily violates 
the local constraint Eq.JSJ. This introduces unphysical 
divergences in correlation functions, which don't have di- 
rect physical meaning. These divergences must cancel for 
all physical properties. According to Read|^, divergent 
goldstone mode phase fluctuations act to restore the true 
symmetry of the system, and (hi) — 0. These phase fluc- 
tuations can be integrated out in the functional integral 
formalism, while the ansatz Eq. (JSJ will hold for the boson 
amplitude. In this sense, this phase transition is similar 
to the Kosterlitz-Thauless phenomenon. In agreement 
with x-ray orthogonality catastrophe [Hj , the local spin 
correlations in the single impurity model are critical, i.e., 
have a power law decay at large times 26] . Fluctua- 
tions in the constraint lead to large deviations from mean 
field at temperatures T ~ Tk- Nevertheless, if the limit 
N — > 00 is taken at finite filling go> the mean field so- 
lution for a single impurity Anderson model approaches 
the Bethe ansatz solution at both high and low tempera- 
ture end |30l ] , with most deviations occurring around the 
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mean field transition temperature T* ~ Tk/IiiN. The 
convergence is excellent at low temperatures, although it 
remains rather poor at T > T* M. This well-known 
problem with the slave boson method has prevented the 
development of the two-fluid microscopic model in the 
past. In particular, the Kondo gas thermodynamic con- 
tribution from almost flat band, which appears already 
in the tree approximation, has always been dismissed as 
small, as an artifact of the above mentioned problem with 
the constraint, Eq.©. 

The mean field Hamiltonian Eq. (JZJ can be diagonalized 
by the canonical transformation, which results in a well 
known 18J spectrum with a gap, 



ei,2(fe) = 2 ( €k + e f ± Ek ) 



E k = J(e k -e f )* + ±(Vaf 



(8) 
(9) 



The new energy spectrum gives rise to the mass enhance- 
ment given by 



m* _ V 2 a z 
m ~ e 2 



1. 



(10) 



which becomes rather dramatic in the Kondo regime. 

It is convenient to introduce re-defined bare-fermion 
Green's functions, 

G™(fc,r) = -(T T c km (T)c{ m (0)) MF , (11) 
G%{k,r) = -(T T f km (r)c{ m (0))MF, (12) 
Gf f (k,r) ee -(T T / fcm (r)/L(0)) MF , (13) 

where the thermal average is taken in a system described 
by Eq.Q and r is imaginary time. Transforming to Mat- 
subara frequencies, it is easy to find: 



G£(fc,« n ) 
GJ f (k,cj n ) 
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- e/)(w„ - 
Va 


- Cfc) 


~V 2 a 2 
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- £fc) 


-V 2 a 2 



(14) 
(15) 
(16) 



The mean field parameters arc determined by minimiz- 
ing the free energy for the Hamiltonian in Eq.Q, which 
diagrammatically leads to vanishing of the sum of the 
"tadpole" graphs |28{. 

N(e f -E )a + VT ]T G%(k,u n ) = (17) 



N(a 2 



qo) 



T E 

k.n.m 



Gf f (k,u n ) = 0, (18) 



where qo is the filling factor, qo — 1/N. For free f- 
electrons above T*, Eq. l|18|l gives: 



n f( e f) = Qo 



l 

N 



(19) 



so that the entropy for the f-levels 

S = -N{q Q In q + (1 - q Q ) In (1 - q )) 



(20) 



Setting the filling q = 1/N gives S ~ In N + 1. Exact 
entropy at high temperatures is not hard to find - one 
has N half-filled bands and In N coming from each iV-fold 
degenerate f- level. Thus, S exac t = In AT, and we make a 
mistake of the order of 1/ In N in the entropy when we do 
mean field. In Coleman approach 26J the 1/N expansion 
takes place at a given qo, and qo is set to 1/N in the end. 
Formally this gives a correction 



Sf, 



S: 



-0.51n(27r7Vg [l - q ] 



(21) 



which is 0(\n N/N), which makes such 1/N expansion 
somewhat better suited for calculating finite temperature 
properties|30j. However, the two-component physics is 
not evident in the Coleman approach. The reason is that 
the Fermi liquid Sommerfeld expansion, which appears 
as an expansion m 1/ In N in the Read and Newns [3 
approach, becomes divergent at T ~ T* in the Coleman 
approach. Hence, the Fermi liquid contribution cannot 
be easily identified. In what follows, we use the approach 
of Read and Newns0| to calculate two-component ther- 
modynamics at the mean field level. 



III. THE MEAN FIELD EQUATIONS. 

The mean field equations EasljT7 )l .l(T5 )l . together with 
the conservation of the total number of particles, 

N tot = T [G/ n /(fc,^„) + G™(fc,c^)] = const, (22) 

k,n,m 

define the order parameter for the second order phase 
transition in the slave boson method. 

Transforming the integration variable in the mean 
field equations Eas H17l) . 118l) to the new energy spectrum 
Eq.©, e = ei,2 (fc), we get|l8|: 



Na 2 



e f - E = p V 
- 1 + PoV 2 



Me) 

(e/-e) 2 



de 



dc. 



(23) 



(24) 



The limits for the integrals are easily calculated from the 
energy spectrum, Eqs||SJ,|J3J: 



-D 



e b = € f 
e c = e f 
e d = D 



D + e f 

(Va) 2 
D-e f 

(Va) 2 

D + €f 

(Vaf_ 



D 



(25) 
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The effective Kondo temperature e/ can be calculated 
from Eq. H23|) . This integral is the Cauchy limit integral 
around the Kondo gap, which is not divergent. In the 
Kondo limit, Eq < 0, \Eq\ ^> V, the zero-temperature 
value of Cf, e/o, is easily calculated from Eo. (|2*3T) : 



De pov' 2 



(26) 



Now, for T ^ 0, but T <C e/, one can use the Sommerfeld 
expansion to get: 

\E \=p V 2 \n 



D 



* T/2 1 

6 ((3e f ) 2 



(27) 



or, equivalently, 



e/ (T)~e /0 (l + — (/3e /0 )- 2 



(28) 



The second mean field equation, Eq.l|24|l. determines 
the order parameter a, or the Kondo gap. The integral 
diverges strongly as (e — e /)~ 2 at e — * £/• To deal with 
this divergence, we can integrate by parts: 



No? 



1 - po 



Yl 

D 



PoD 



D 



+ n f \€ f + 



r 



PoV 2 



n'f(e)de 



(29) 



The latter integral is now convergent in Kondo gap re- 
gion, so it can be taken using Sommerfeld expansion, 



PoV- 
To solve for T 

which gives 



n' f {e)de Po V 2 



we take the limit a 2 - 
q = 2p Q Dn f (e f ), 



1 



TT 2 T 2 

3ej 



(30) 



in Eq.^: 



In [2 P() DN] 



(31) 



(32) 



The zero-temperature value of the order parameter a, 
a(T = 0) = ao, can also be easily found from Ea. H29|) : 



2 



1 



PoV 2 



(33) 



TV e /0 

The transcendental equation Ea. H29|l . which gives the 
temperature dependence of the order parameter a, is 
considerably simplified if we only calculate the lowest 
terms in 1/lnTV. Indeed, since T < T* ~ e//lniV, the 
asymptotic Sommerfeld expansion becomes an expansion 
in 1/lniV. To keep the lowest order in 1/lnTV in the 
Sommerfeld expansion, we take 



e f (T < T*) ~ e /0 



(34) 



in Ea. H32|l . One can simplify the gap equation Eq. l|29|l in 
the Kondo regime, where 



A P° V2a2 ^ 
A fc = D < e /- 



(35) 



Assuming a square band 2poD — 1, following Newns and 
Read fl8l|. keeping the leading terms in 0(1/ IniV), 



a 2 (T) 



l-Nn f (e f ) 



(36) 



The third equation, Eq. I122H , determines the shift of the 
chemical potential in the slave boson phase transition. 
In the Kondo limit the number of electrons per /-level 
approximately stays at 1, and no significant shift of the 
chemical potential occurs. 



IV. THERMODYNAMICS OF THE KONDO 
LATTICE - THE TWO-COMPONENT PICTURE. 

Thermodynamic properties can easily be calculated 
from the mean field free energy for the Hamiltonian 
Eq.Q. From Luttinger functional, we can write for the 
Grand free energy density without the small free conduc- 
tion electron contribution: 

fi = -P = ^2n f + Sn + ^2N{a 2 -q )(ef-E ), (37) 



where flf is the grand free energy density for free non- 
hybridized f-centers, Sft is the effect of hybridization 
(mixed contribution). Since flf corresponds to a free 
fermion at renormalized energy e / with respect to chem- 
ical potential, 



Q f = -NT\n(l + e 



(38) 



The third term in Eq. l|37|l comes from the energy shift 
due to a slave boson average. The interaction energy <5f2 
can be expressed through the mixed Green's function, 



Sn = 2V 



(39) 



fc,n,m 



We note that the RHS of this equation corresponds to 
the Green's function expression in Eq. l|18|l . We can 
use Eq. l|18|) to replace the expression under the integral. 
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Since from Eq. itTBl) 

VT G^(iu n ,k) = -c f N(e f - Eo)a, (40) 

k.n,m 

Sn = -2c f N [ (e f -E )ada = -Nc f a 2 (e f - E ), (41) 
Jo 

where c/ is the density of f-levels. As a result, we easily 
obtain: 

fl = -NcfTln (1 + e-^f) - Nc f q e f (42) 
I 



Comparing this with expression Ea. ()^6jl for a 2 , and keep- 
ing the terms of lowest order in 0(1/ In TV) (the first term 
is then small compared with the last term), we can write: 

2 _ 2 2 y2 2 

S = Pe f c f a -^+\N Po V -^T, (44) 

or 

S~~/ HFL Tf(T)+Pe f c f (l-f(T)), (45) 
I 



Houghton, Read, and Won|33( have found that a similar 
expression for the Grand free energy for a single- impurity 
model holds to two loops. Thermodynamics is obtained, 
as usual, by differentiating the Grand Free energy with 
respect to temperature: 



(43) 

I 

where 



f(T) = , (46) 



and jhfl — ^r-Npo^. The heat capacity can be 
found by differentiating the entropy S: 



S= -^T = NCf ln (1 + 6 ~^ e/) + (1 ~ Nn f( € f» + fcfNc f n f ( ef ) 



dS e 2 a 2 , - a 2 vr 2 V 2 a 2 



where we have neglected all terms smaller in T/ej ~ 
1/lnTV in the first part, although all of them can be 
written out in the same fashion: 



r-T dS -r r a °~' 
C - T dT - ClCf ^ 



Cfl- 



(48) 



"0 "0 

Here Cfl is the usual heavy Fermi liquid contribution, 
C FL = IhflT, (49) 

while 

Ci - ^4 , (50) 

a factorized contribution. We note that unlike Ref . [l2l| . 
the dominant contribution in 1/ In N comes from differ- 
entiating the order parameter a 2 (T), 

d(a 2 ) 



dT 



J^2 



{a 2 -a 2 ). 



(51) 



In magnetic field we can write, neglecting free conduc- 
tion electron contribution: 



Q = -Tc f J2 ^ (l + e-^/+ ftm )J - c f e f (52) 

m 

Now, N in MF equations is replaced by sum over different 
ef(h,T) + hm. Taking a partial derivative with respect 
to h — g^sB, we find the magnetic moment: 



M 



9HB 
For small h 



E ( + m ) n f( e f + hm ) +c f^h (53) 



dc 



= c f^O--Nn f (e f ))-hc f J £m l n' / (e f ) (54) 

This obviously separates into two contributions. After 
some simple algebra, we find that the bulk spin suscep- 
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Kondo gas 




FIG. 1: The heavy liquid component comes from the renor- 
malized band contribution, while the gas component appears 
due to thermal population of almost localized part of the spec- 
trum. 



corresponds to the fraction of the Kondo liquid, Ea. (|46|) . 
The linear behavior of f(T) for T < T* , which is dif- 
ferent from that given by Eg. 1361) . was found by fitting 

Order parameter for N = 6 (Ce) 




T/T* 



tibility also has two separate contributions: 



a 2 . (9VB) 2 j(j + l)c/ 
X - Xhfl — H 



3T 



1 - 



where 



NiWBfjjj + 1) rn* 

Xhfl 5 Po — 

3 m 



(55) 



(56) 



and m*/m is the mass enhancement at T = 0, j = ~ . 
The second term on the rhs corresponds to the typical 
Curie impurity contribution. A correct logarithmic cor- 
rection to it appears in the next order in l/jVj^ 

The origin of the two components in thermodynamics 
is shown schematically on Figd Indeed, the new energy 
spectrum for free fermions gives a singularity in density 
of states, 



= Po 



1 



(Va) 



(57) 



This singularity contributes to thermodynamic proper- 
ties. Since the mean field equations are derivatives of the 
free energy, the singularities that are encountered are of 
the same type as those for ef and a 2 . Therefore, the final 
results can be expressed in terms of a 2 (T). 



V. SUMMARY 

Comparing the mean field results with experimental 
phenomenology of Ref. |T^ . we find many similari- 
ties. For example, thermodynamics for the Kondo lattice 
in the slave boson approach clearly factorizes into the 
Kondo impurity and heavy electron liquid contributions. 
Similar to situation in superconductivity this factoriza- 
tion is a direct result of the appearance of anomalous 
mixed averages of the type (cr /), which produces the en- 
ergy gap, or the mean field order parameter a 2 , which 



FIG. 2: The temperature dependence f(T) = 
(N=6) 



for Ce 



experiments in the 1-1-5 family in Ref. [T^| . In case of 
Ce (N — 6) we find that numerically Ea. H36l) gives an 
order parameter which is almost linear (See Fig[2J) The 
spin susceptibility x(T) follows the experimental behav- 
ior, and the behavior of thermodynamic functions factor- 
izes into two components as well. Clearly this is a very 
rough approximation, and many effects observed exper- 
imentally in Ref.|l2|| are not captured. For example, at 
zero temperature we find that f(T = 0) = 1, and the 
Kondo impurity contribution is absent, unlike the limit 
f(T = 0) ~ 0.9 found in Ref. [HI. This discrepancy 
could, in principle, be explained by disorder, which would 
produce a finite population of the DOS peak at zero tem- 
perature. We also note that the above model is consider- 
ably oversimplified. For example, the fc-dependence of V 
was neglected. In a realistic case V(k) may have zeroes 
on the Fermi surface, which will lead to a fc-dependent 
energy gap with a nodal structure, similar to that ob- 
served in unconventional superconductors. 

We also find that thermodynamics (C, S, etc.) does 
not completely follow the NPF phenomenology. While 
the separation in two components is always present, we 
found that it is necessary to differentiate f(T) to get the 
correct answer for some quantities. It is also, perhaps, 
necessary to get to the next approximation, in order to 
find the log(T) dependence of effective mass, observed 
in Ref.|l2l|. Calculation of transport properties requires 
going beyond the mean field approximation |28| . which, 
in principle, could also explain f(T = 0) ^ 1 without 
introducing disorder. Nevertheless, this approach can be 
regarded as the first approximation to understanding ex- 
perimental results. 

I would like to thank L.P. Gor'kov for many discussions 
and insights. This work was supported by TAML at the 
University of Tennessee. 
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